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Abstract 

Bacteria rely on two known DNA-level defenses against their bacteriophage predators: restriction-modification and 
Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-CRISPR-associated (Cas) systems. Certain phages have 
evolved countermeasures that are known to block endonucleases. For example, phage T4 not only adds hydroxymethyl 
groups to all of its cytosines, but also glucosylates them, a strategy that defeats almost all restriction enzymes. We sought to 
determine whether these DNA modifications can similarly impede CRISPR-based defenses. In a bioinformatics search, we 
found naturally occurring CRISPR spacers that potentially target phages known to modify their DNA. Experimentally, we 
show that the Cas9 nuclease from the Type II CRISPR system of Streptococcus pyogenes can overcome a variety of DNA 
modifications in Escherichia coli. The levels of Cas9-mediated phage resistance to bacteriophage T4 and the mutant phage 
T4 gt, which contains hydroxymethylated but not glucosylated cytosines, were comparable to phages with unmodified 
cytosines, T7 and the T4-like phage RB49. Our results demonstrate that Cas9 is not impeded by N'^-methyladenine, 5- 
methylcytosine, 5-hydroxymethylated cytosine, or glucosylated 5-hydroxymethylated cytosine. 
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Introduction 

Bacteria utilize an assortment of anti-phage defense mecha- 
nisms, including two that act at the nucleic acid level: restriction- 
modification and Clustered Regularly Interspaced Short Palin- 
dromic Repeats (CRISPR)-CRISPR-associated (Cas) systems. 
Some bacteriophages have developed extensive modifications to 
their DNA that enable them to evade host restriction endonucle- 
ases. For example, phage T4 replaces each cytosine with 
hydroxymethylated cytosine (hmC), then glucosylates the hydro- 
xymethyl group to form glucosylated hmC (ghmC) [1]. The bound 
glucose shelters the phage genome from the host's modified 
cytosine restriction systems, McrA, McrBC, and Mrr, which 
recognize methylcytosines and hmCs but not ghmCs [2]. 

CRISPR-Cas systems also function as endonucleases, though 
unlike restriction enzymes, their recognition sites are programma- 
ble by CRISPR RNAs (crRNAs) [3]. As an adaptive immune 
system, CRISPR-Cas components incorporate fragments of DNA 
from invading viruses or plasmids into arrays composed of spacers 
interspersed with repeats on the genome [4,5]. In Type II 
CRISPR systems, transcribed arrays are processed into crRNAs 
that form a complex with the RNA-guided Cas9 nuclease and a 
trans-activating RNA (tracrRNA) [6]. The crRNA guides the 
complex to double-stranded DNA "protospacer" sequences that 
match the sequence of the spacer and are flanked by a 
"protospacer adjacent motif (PAM) unique to the CRISPR 



system [7]. If spacer-protospacer base-pairing is a close match, 
Cas9 cuts both strands of DNA, often eliminating the plasmid or 
phage. We sought to determine whether various DNA modifica- 
tions known to block restriction systems can similarly impede 
CRISPR-Cas defenses. 

Materials and Methods 

Bioinformatics search 

We derived a list of 1749 unique spacers from several sources: 
49 E. coli strains with CRISPR structures in the CRISPRdb 
database (http://crispr.u-psud.fr/crispr/, [8]), 72 strains in the 
ECOR collection [9], 263 strains isolated from humans or animals 
in various regions of France [10], and 194 Shiga toxin-producing 
E. coli (STEC) strains [1 1] . CRISPR array sequences were 
processed in CRISPRfinder (http://crispr.u-psud.fr/Server/, 
[12]) to extract spacer sequences. 

We performed BLASTn searches (http://blast.ncbi.nlm.nih. 
gov/, [13]) with a word size of seven optimized for short sequences 
and an E-value of less than 0.1, which corresponded to roughly at 
least 14 matched nucleotides in the T2/T4/T6 genomes search 
and at least 17 matched nucleotides in the all T4-like genomes 
search. We screened hits by first looking for a concentration of 
exact nucleotide matches at the 5' end, which would be consistent 
with a seven-nucleotide "seed" region that does not tolerate 
mismatches [14]. Outside the seed sequence, at least five 
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mismatches are tolerated [14], though the upper limit of tolerable 
mismatches has not been characterized in the E. colt CRISPR 
system. We then checked for a properly oriented E. coli Type I-E 
CRISPR PAM such as AAG, ATG, AGG, and GAG in the 

targeted sequence. 

Bacterial strains and plasmid construction 

In addition to wild-type E. coli K-12 MG1655 and E. coli B, we 
used methyltransferase-deficient (dam~ / dcm~) E. coli K-12 
(ER2925, New England Biolabs, Ipswich, MA) and restriction- 
deficient [mcrA' mcrBCr mrr~ hsdRT) E. coli K-12 (ER1821, New 
England Biolabs). E. coli were grown at 37°C in LB broth and 
supplemented with antibiotics as needed at final concentrations of 
100 |Xg/mL spectinomycin, ,30 |Xg/mL chloramphenicol, 300 Hg/ 
mL erythromycin, and 100 |tg/mL carbenicillin. 

Cells expressing SpCas9 were constructed by transforming in 
DS-SPcas (Addgene plasmid 48645, [15]), which encodes SpCas9 
and its cognate tracrRNA on a backbone with a cloDFlS origin of 
replication and aadA gene. In the dam/dcm methylation studies, 
we assembled a compatible protospacer plasmid encoding all five 
of the target sequences with their PAMs; we placed the control, 
daml, and dcml sequences after a pBR322 origin of replication, 
and the dam2 and (lcm2 sequences after a bla gene. In the T7 
infection assays, the spacer was expressed on DS-SPcas such that 
there was no separate spacer plasmid. In all other experiments, we 
maintained the designed spacer on a separate plasmid (based on 
PM-SPITB, Addgene plasmid 48650, [15]) that expressed one 
spacer followed by the SpCas9 repeat on a backbone with a pi 5a 
origin of replication and cat gene. When a dilferent resistance 
marker was needed, we switched cat with EryR. 

Bacteriophage strains and propagation 

Phage T7 stock was propagated in E. coli K-12 MG1655 and 
RB49 stock (obtained from H. M. Krisch) propagated in E. coli B. 
Wild-type T4 stock was propagated in E. coli K-12 MG1655. 
Phage T4 gt (a gift from New England Biolabs) is T4 (».-gt57 p-gtl4, 
which does not have functional a- and P-glucosyltransferases [16]. 
Because the E. coli restriction system recognizes and cleaves hmC, 
preventing T4 gt from plaquing efficiendy, we conducted all 
experiments involving this phage in the restriction-deficient E. coli 
K-12 host ER1821. 

In phage stock preparation, an overnight bacterial host culture 
was diluted 1:100 in LB, inoculated with phage, and grown for 
2.5-5 hours (during which the turbidity of cultures rose and then 
fell due to lysis). The lysates were spun down at 8000 x g for 5 
minutes at 4°C to remove cell debris. The supernatant was filtered 
through a 0.45 nM membrane and stored at 4°C. 

Transformation assays 

We prepared protospacer and spacer plasmids from a dam"""/ 
dcm""" strain, NEB Turbo (New England Biolab.s), and performed 
transformation assays using E. coli K-12 MG1655 bacteria 
containing the protospacer plasmid and DS-SPcas. After trans- 
forming equimolar amounts of each spacer plasmid and selecting 
for all three plasmids (DS-SPcas, protospacer, and spacer), we 
quantified the number of transformants relative to a transformed 
spacer plasmid that did not target the protospacer plasmid. We 
also reversed the transformation order for one set of experiments; 
that is, we transformed the protospacer plasmid into E. coli already 
carrying DS-SPcas and each spacer plasmid. We observed 
comparable numbers of transformants regardless of order. We 
repeated the same transformations in methyltransferase-deficient 
E. coli K-12 using equimolar unmethylated protospacer and spacer 
plasmids, which were prepared from E. coli K-12 damT IdcmT . 



Again, for one set of experiments, we reversed the transformation 
order and noted similar numbers of transformants. 

Plaque assays and efficiency-of-plating calculations 

To characterize the level of phage resistance conferred by Cas9, 
we infected normalized densities of protected E. coli with equal 
titers of phages and counted the number of plaques. Equal cell 
densities were obtained by diluting an overnight culture and 
normalizing to an ODgooj,^ of 0.3 after several hours of growth. 
We added 2 |J,L of phage to 120 |J,L of cells, mixed them 
thoroughly in 1 mL of 0.6% top agar with appropriate; antibiotics 
within 20 minutes, and poured the mixture onto 3 mL of 1.5% 
sohd agar. Independent experiments were performed with 
dilferent phage dilutions. To calculate an efficiency of plating 
(EOP), we divided the phage titer from plating the phage on a 
protected strain by the phage titer from plating the phage on a 
susceptible wild-type strain. 

Results 

Natural spacers target phages with modified DNA 

We began by attempting to discover naturally acquired spacers 
in bacteria that target phages known to contain modified DNA. 
Only a handful of phage families have been identified with 
completely modified DNA, including Bacillus subtilis phage PBS2, 
Synechococcus elongates phage S2L, and Eschemhia coli phage T4 [17]. 
Since CRISPR-Cas systems and phages of E. coli have been better 
studied than those of the other bacterial hosts, we focused on 1749 
unique E. coli spacers in available array sequences from the ECOR 
collection, Shiga toxin-producing E. coli (STEC), and other 
databases. 

Upon searching for candidate protospacers in phages T2, T4, 
and T6, all of which contain ghmC DNA [1], we found one hit 
that matched 25 of 32 nucleotides in T2's gene 38, although this 
spacer was only found in one human-associated E. coli (Figure 1 A). 
In an expanded search including T4-like phages, we identified 
another hit with 29 nucleotides matching phage CC3rs gp35 
(Figure IB). CC31 is the only known non-T-even type phage with 
predicted glucosyltransferase genes [18], which are required for 
generating ghmC from hmC. This spacer was found in many 
different E. coli isolates. 

The potential presence of natural spacers targeting phage with 
modified DNA suggests that CRISPR-Cas systems may overcome 
this form of phage defense. To test this hypothesis, we explored the 
extent to which the Type II-A Streptococcus pyogenes Cas9 (SpCas9), 
the most commonly used CRISPR-Cas system for genome 
engineering, is able to cleave various forms of modified DNA. 

Cas9 cuts N^-methyladenine and 5-methylcytosine in 

E coli 

DNA adenine methyltransferase (dam) methylates the adenine 
in 5'-GATC-3', while DNA cytosine methyltransferase (dcm) 
methylates the internal cytosine in 5'-CCTGG-3' and 5'- 
CCAGG-3' in E. coli. We designed target sequences containing 
one to two dam or dcm sites as well as a control target sequence 
with no methylation sites (Figure 2A). We prepared spacer and 
protospa( (;r plasmids from a dam''^/dcm^ strain and selcct(;d for 
the coexistence of each spacer and its targeted protospacer in 
transformation assays using dam^/dcm^ cells expressing SpCas9. 
AU targeted sequences yielded 10^ to 10"* fewer transformants than 
the non-targeted control regardless of whether they contained dam 
or dcm methylation sites (Figure 2B). We observed similar values 
in methyltransferase-deficient (dam Idem ) E. coli K-12, in which 
all plasmids were prepared from a damT /dcm~ strain and were 
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spacer S641 

5 ...CCGACGCATTCGAACTGGGCGACATAGAGCCAGGCGAG...3' 

3'...GGCTGCGTAAGCTTGACCCGCTGTATCTCGGTCCGCTC...5' 



E. coli CRISPR array 



B 



5' CCGACGCAUUCGAACUGGGCGACAUAGAGCCAGGCGAG 3' CrRNA 

IIIIIIMIIII II I I IMMI III 

3'...TTGTGCGTAAGCTTGGCCTGTTATATCTCAGTCTCGGA.„5' 
5 ...AAGACGCATTCGAACCGGACAATATAGAGTCAGAGCCT...3' 

PAM protospacer 



Phage 12 



spacer S134 

5 ...CCGATTTCATCAAAGCATTAAGGGATGGAATAAAGGAG...3' 

3'...GGCTAAAGTAGTTTCGTAATTCCCTACCTTATTTCCTC,..5' 



E. coli CRISPR array 



5' CCGAUUUCAUCAAAGCAU UAAGGGAUGGAAUAAAGGAG 3' CrRNA 

I I I I I I I I I I I I I I I Mill I llllllll 

3'...GTCTAAAGTAGTTTCGTAGATGTCTTCCCGAACTTATTTCGGA...5' 
5'...GAGATTTCATCAAAGCATCTACAGAAGGGCTTGAATAAAGCCT,.,3' 

PAM protospacer 



Phage CC31 



Figure 1. Native £ co// spacers target phage witli modified DNA. In a BLASTn search, 1749 unique spacers from sequenced £ coli CRISPR 
arrays were queried against T4-like phage genomes. (A) Spacer S641 matches 25 of 32 nucleotides in phage T2. The putative protospacer has a 
permissible £ coli CRISPR PAM AAG and the matching nucleotides are concentrated at the 5' end as a seed sequence. The spacer originated from the 
CRISPR1 locus of £ coli strain 579, a human-associated isolate from France. (B) Spacer SI 34 matches 29 of 32 nucleotides in phage CC31 . While the 
protospacer in phage CC31 has five nucleotides inserted in the center of the sequence, there are 15 exactly matched nucleotides at the 5' end in 
addition to 14 matched nucleotides after the insertion. The PAM GAG and strongly matched seed region suggest it is a plausible £ coli CRISPR target. 
This spacer was found in several strains, including £ coli C str. ATCC 8739, ECOR strains 17 through 21, one farm pig and two human fecal samples in 
France, duck and cattle fecal samples in Australia [33], and enterotoxigenic £ coli (ETEC) strain UMNK88. The spacer and matching protospacer are in 
blue, the transcribed CRISPR RNA (crRNA) in bold black, and PAM sequence in red. 
doi:10.1371/journal.pone.0098811.g001 



thus unmethylated. Overall, we detected no difference in Cas9 
activity on adenine-methylated, cytosine-methylated, and un- 
methylated target sequences. These results are consistent with 
reports showing adenine methylation does not affect CRISPR- 



mediated phage resistance in Streptococcus thermophilus [19] and 
cytosine methylation does not affect SpCas9 activity on sequences 
with CpG sites in human cells [20]. 
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Figure 2. Cas9 cuts metKiylated cytosines and adenosines in £ coli. (A) Synthetic targets were designed to contain one to two dam (orange) 
or dcm (blue) sites. A control unmethylated sequence (+) was included. The PAM sequence NGG for SpCas9 recognition is underlined. (B) In serial 
transformations, we selected for the coexistence of DS-SPcas, the protospacer plasmid, and each spacer plasmid. The number of transformants was 
divided by the number of colonies resulting from a control transformation using a spacer plasmid (-) that did not target the protospacer plasmid. This 
relative number of transformants is plotted for £ coli K-1 2 and £ coli K-1 2 dam^ /dcm" from three independent experiments. Lines represent the 
median. 

doi:1 0.1 371 /journal.pone.009881 1 .g002 
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Cas9 provides resistance against phages T7 and T4-like 
RB49 with unmodified DNA 

We next tested the ability of SpCas9 to provide resistance to 
lytic phages without DNA modifications by constructing spacers 
against phages T7 and RB49, neither of which contains modified 
DNA. RB49 is a T4-lik:e phage that is missing hydroxymethylase 
and P-glucosyltransferase, which are required for modifying 
cytosine to hmC and hmC to P-ghmC, respectively [21]. We 
designed four spacers: two targeting the gene encoding the 
primase/helicase enzyme of T7 (Figure 3A) and two targeting the 
gene encoding the major capsid protein of RB49 (gp23), which is 
one of the most conserved regions across T-even phages [21] 
(Figure 3C). We transformed each spacer-encoding plasmid into 
SpCas9-expressing E. coli K-12 MG1655 and E. coli B to create 
strains protected from T7 and RB49 infection. We challenged 
these strains with phage to calculate an efficiency of plating (EOF) 
compared to unprotected strains; representative plaque plates are 
included (Figures 3A and 3C). 

In E. coli B, T7 had an EOF of 10 on cells expressing spacer 1 
or 2 relative to cells without spacers (Figure 3B). In E. coli K-12, 
spacer 1 reduced sensitivity to T7 infection by four orders of 
magnitude, though spacer 2 only lowered sensitivity by one order 
of magnitude for unknown reasons. RB49 had an EOF of 1 0 ''on 



E. coli B with spacer 1 or 2, and an EOF of 10 ' on E. coli K-12 
with spacer 1 or 2 (Figure 3D). The decreased plaquing efficiencies 
of T7 and RB49 on protected strains reflect Cas9 activity against 
invading unmodified phage DNA. 

Cas9 provides resistance against mutant phage T4 with 
hmC DNA and wild-type T4 with ghmC DNA 

Having established that Cas9 can confer resistance against non- 
modified phage, we proceeded to challenge it with T4 phage 
containing either hmC or ghmC DNA. During replication, wild- 
type T4 synthesizes hmC, which contains a hydroxymethyl group 
attached to the C5 position of cytosine, by using hydroxymethy- 
lated dCTF serially converted from dCTF [22]. Then phage- 
encoded glucosyltransferases add a glucose group to the hydro- 
xymethyl group in ot- or P-configuration [23] (Figure 4A). To 
investigate Cas9 activity against T4 without glucosylated DNA, we 
included mutant phage "T4 gt", which has hmC rather than 
ghmC due to non-functional glucosyltransferases [16]. By using 
restriction enzymes with varying sensitivity to modified cytosines 
(according to REBASE, http://rebase.neb.com/), we confirmed 
that our stocks of phage T4 had ghmC, phage T4 gt had hmC, 
and phage RB49 did not have ghmC or hmC (Figure SI). 
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Figure 3. Cas9 reduces £ co// susceptibility to phages T7 and RB49. (A) Spacers against T7 were targeted against the primase/helicase gene 
(gene 4A and 4B). The PAM is underlined in the sequence and shown as a black box in the diagram showing the orientation and location of the 
protospacer (white box) on the gene. In a representative T7 plaque assay of protected and unprotected strains, there is substantial lysis on wild-type 
(wt) £ coli K-12, visible plaquing on cells with spacer 2 (sp 2), and no plaques on cells with spacer 1 (sp 1). (B) The efficiency of plating of T7 was 
calculated for each protected strain relative to the unprotected wild-type strain. Independent replicates of £ coU B (n = 4, 3, 3) and £ coli K-12 (n = 5, 5, 
7) are plotted. Lines represent the median. (C) Spacers against RB49 were constructed against the major capsid protein (gp23). In a typical RB49 
plaque assay, there is notable lysis on wild-type £ coli B, some plaques on cells with spacer 1, and a few plaques on cells protected with spacer 2. (D) 
The efficiency of plating of RB49 was quantified for each protected strain relative to the unprotected wild-type strain. Shown are independent 
replicates of £ coli B (n = 5, 3, 3) and £ coli K-12 (n = 3, 3, 3). Lines represent the median. 
doi:1 0.1 371/journal.pone.009881 1 .g003 
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Figure 4. Cas9 reduces £ co//susceptibility to phages T4 and T4 gt. (A) The structures of cytosine and modified cytosines are sinown. T4 gt 
has 100% hydroxymethylated cytosines (hmCs). T4 has 100% glucosyl-hydroxymethylated cytosines (ghmCs), specifically 70% at- and 30% |3-ghmCs. 
The ghmC structure shown is in the p-configuration. (B) Spacers against T4 were also designed against the major capsid protein (gp23), which is 
homologous to that of RB49. For comparison, the RB49 protospacers are aligned below in italics, where dots indicate identical nucleotides. In the T4 
sequences, the PAM is underlined. The PAM (black box) and protospacer (white box) are represented on the gene. (C) In a typical plaque assay with 
T4 gt (left plate), there was complete lysis on wild-type (wt) restriction-less (r-l) £ coli K-12 and few plaques on cells with spacers 1, 2, or 3 (sp 1, sp 2, 
or sp 3). In an assay with T4 (right plate), there was complete lysis on wild-type £ coli K-12 MG1655, numerous plaques on cells with spacer 1 or 3, and 
about a dozen on spacer 2. (D) The efficiency of plating of T4 and T4 gt was quantified for each protected strain relative to the unprotected wild-type 
strain. Independent replicates of restriction-less £ coli K-12 (n = 5, 3, 3, 5), £ coli K-12 (n = 4, 4, 5, 6), and £ coli B (n = 5, 3, 3, 3) are plotted. Lines 
represent the median. 
doi:1 0.1 371 /journal.pone.009881 1 .g004 



Since T4 gp23 is liomologous to gp23 from RB49, we modified 
our two spacers against RB49 to matcli the sequences of T4, and 
also designed an additional spacer (Figure 4B). We tested these 
spacers using efficiency-of-plating experiments as before; repre- 
sentative plaque plates are shown (Figure 4C). Assays involving T4 
gt used restriction-less E. coli K- 1 2 because wUd-type K-12 restricts 
hmC DNA; the EOF of T4 gt on E. coli K-12 MG1655 is 10"* 
compared to T4 on MG1655 (Figure S2). In restriction-less E. coli 
K-12, T4 gt exhibited an EOF of 10"*^ to 10"'' on ceUs carrying 
any one of the three spacers (Figure 4D). Wild-type T4 displayed 
an EOF of lO""^ on E. coli K-12 MG1655 with spacers 1 or 2, and 
an EOF of 10~ ^ on cells expressing spacer 3. On E. coli B with any 
three spacers, T4 had an EOF of 10 to 10 *. As the dilference 
in EOF values for both T4 gt and wild-type T4 phages were 
comparable to those of the non-modified T4-like phage RB49, our 
results demonstrate that SpCas9 is not impeded by hydroxy- 
methylation or glucosyl-hydroxymethylation of phage DNA. 

Discussion 

Our discovery that S. pyogenes Cas9 is insensitive to methylation, 
hydroxymethylation, and glucosyl-hydroxymethylation renders it 
unique among current genome-targeting technologies, as both 
zinc-fingers (ZFs) and transcription activator-like (TAL) effectors 
can be engineered to discriminate 5-methylcytosrne from cytosine 



[24,25]. This difference may be useful for biotechiiological 
apphcations. 

In our bioinformatics search for candidate natural spacers, we 
were only able to identify two possible sequences against T4-like 
phages. This type of bioinformatics search is hampered by the 
currently limited knowledge of specificity and tolerability of 
mutations in both the acquisition and interference stages of 
CRISFR systems. While this paper was under review, Fineran et 
al. published a report exploring the robustness of the E. coli 
CRISFR system, in which degenerate target regions with up to 13 
mutations in the protospacer and FAM can promote "priming," a 
positive-feedback mechanism to incorporate new spacers based on 
mutated or outdated spacers [26]. This suggests more lenient 
bioinformatics searches would be allowable. Furthermore, our 
search is limited by available sequences of E. coli and phages 
known to modify their DNA, as well as the possibility that these 
isolates do not encounter T4-like phages in their environments. 
Future searches may provide additional evidence of CRISFR- 
based immunity to DNA-modifying phages. 

Interestingly, we observed that different spacers conferred 
differing levels of resistance against phage infection. Since 
mutations in the protospacer or FAM can allow phage to escape 
[27,28], we sequenced Cas9-targeted regions of plaques that 
appeared on protected strains. Indeed, T4 and T7 plaques on 
protected E. coli had mutated one nucleotide in the FAM, or one 
to two nucleotides in the protospacer (Table SI). Less effective 
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spacers may be targeting sequences that are more readily mutated, 
though we cannot rule out the non-mutually exclusive possibility 
that Cas9 acts more slowly on certain sequences and thus allows 
phagp-induced lysis to outpace Cas9-enabled protection. In S. 
thermophilus CRISPRl and CRISPR3 systems, the uncut phage 
genome can stiU be observed in bacteriophage-insensitive mutants 
[29,30]. Further investigation of how some but not all phage DNA 
molecules escape Cas9 cutting during phage infection is needed. 

While phages may inactivate CRISPR proteins [31] or encode 
their own CRISPR-Cas systems [32], we have demonstrated that 
DNA modifications that normally circumvent bacterial restriction 
systems do not impede Type 11 CRISPR systems. Our findings 
may help explain why DNA modifications remain uncommon 
among bacteriophages characterized to date whereas nearly half of 
bacteria have CRISPR structures [8]. 

Supporting Information 

Figure SI Restriction digest of phages. Phage DNA was 
extracted by using the Qiagen Blood and Tissue Kit on 200 (iL of 
phage stock. 10 or 20 U of each enzyme (1 |j.L) was added to 5 \xL 
of lOX CutSmart Buffer (NEB) in a 50 |tL reaction volume 

containing approximately 1 00 ng of phage RB49 or T4 DNA, or 
800 ng of T4 gt DNA. The reactions were incubated at 37°C for 
4 hours before visualizing on a 1 % agarose gel stained with SYBR 
Gold. As expected, Dral cuts all RB49, T4 gt, and T4; Hpall and 
Nhel are sensitive to methylated cytosines and only cut RB49; and 
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Xbal has 50% activity on hmC and partially cuts T4 gt. Blue text 
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described in the main text. 
(TIF) 
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are underlined. Mutations are in bold text and double-underlined. 
(DOCX) 

Author Contributions 

Conceived and designed the experiments: SJY KME. Performed the 
experiments: SJY. Analyzed tlie data: SJY KME CMC. Contributed to the 
writing of the manuscript: SJY KME CMC. 



15. Esvelt KM, Mali P, Brafl'JL, Moosburner M, Yaung SJ, et al. (201 3) Orthogonal 
Cas9 proteins for RNA-guided gene regulation and editing. Nat Methods 10: 
1116-1121. doi:10.1038/nmedi.2681. 

16. Georgopoulos CP (1967) Isolation and preliminary characterization of T4 
mutants with iiongliicosylated DNA. Biochem Biophys Res Gommun 28: 179— 
184. doi:http://dx.doi.o'rg/10.1016/0006-291X(67)90426-3. 

17. Warren RAJ (1980) Modified bases in bacteriophage DNAs. Annu Rev 
Microbiol 34: 137-158. doi: 10.1 146/annurev.mi.34.100180.001033. 

18. Petrov VM, Ratnayaka S, Nolan JM, Miller ES, KaramJD (2010) Genomes of 
the T4-related bacteriophages as windows on microbial genome evolution. 
VirolJ 7: 292. doi:10.1186/1743-422X-7-292. 

19. Dupuis M-E, Vfflion M, Magadan AH, Moineau S (2013) CRISPR-Gas and 
restriction-modification systems are compatible and increase phage resistance. 
Nat Commun 4: 2087. doi:10.1038/ncomms3087. 

20. Hsu PD, Scott DA, VVeinstein JA, Ran FA, Konermann S, et al. (2013) DNA 
targeting specificity of RNA-guided Cas9 nucleases. Nat Biotechnol 31: 827- 
832. doi:10.1038/nbt.2647. 

21. Monod G, Repoila F, Kutateladze M, Tetart F, Krisch HM (1997) The genome 
of the pseudo T-even bacteriophages, a diverse group that resembles T4. J Mol 
Biol 267: 237-249. doi:10.1006/jmbi.l996.0867. 

22. Snyder L, Gold L, Kutter E (1976) A gene of bacteriophage T4 whose product 
prevents true late transcription on cytosine-eontainiiig T4 DNA. Proc Natl Acad 
Sci USA 73: 3098-3102. 

23. Kornberg S, Zimmerman S, Komberg A (1961) Glucosylation of deoxyribo- 
nucleic acid by enzymes from bacteriophage-infected Escherichia coli. J Biol 
Chem 236: 1487-1493. 

24. Choo Y (1998) Recognition of DMA mcthylation by zinc fingers. Nat Struct Biol 
5: 264-265. doi:10.1038/nsb0498-264. 

25. Valton J, Dupuy A, Daboussi F, Thomas S, Marcchal A, et al. (2012) 
Overcoming transcription activator-like effector (TALE) DNA binding domain 
sensitivity to cytosine methylation. J Biol Chem 287: 38427-38432. 
doi: 10.1 074/jbc.C 1 1 2.408864. 

26. Fineran PC, Gerritzen MJH, Suarcz-Dicz M, Kunne T, Boekhorst J, et al. 
(2014) Degenerate target sites mediate rapid primed CRISPR adaptation. Proc 
Natl Acad Sci U S A. Available; http://w\\'w.ncbi.nlm. nih.gov/pubmed/ 
24711427. Accessed 9 April 2014. 

27. Deveau H, Barrangou R, Garneau JE, Labonte J, Fremaux C, et al. (2008) 
Phage response to GRISPR-encoded resistance in Streptococcus thermophilus. 
J Bacteriol 190: 1390-1400. doi:10.1128/JB.01412-07. 

28. Levin BR, Moineau S, Bushman M, Barrangou R (2013) The population and 
evolutionary dynamics of phage and bacteria with CRISPR-mediated immunity. 
PLoS Genet 9: el003312. doi:10.1371/journal.pgen.l003312. 

29. Magadan AH, Dupuis M-E, Villion M, Moineau S (2012) Cleavage of phage 
DNA by the Streptococcus thermophilus CRISPR3-Cas system. PLoS One 7: 
e409 1 3. doi: 10. 1 37 l/joumal.pone.00409 1 3. 



PLOS ONE I www.plosone.org 



6 



June 2014 I Volume 9 | Issue 6 | e98811 



DNA Modifications Do Not Impede CRISPR/Cas9-Mediated Phage Resistance 



Gamcau JE, Dupuis M-E, ViUion M, Romero DA, Barrangou R, ct al. (2010) 
The CRISPR/Cas bacterial immune system cleaves bacteriophage and plasmid 
DNA. Nature 468: 67-71. doi:10.1038/nature09523. 

Bondy-Denomy J, Pawluk A, Maxwell KL, Davidson AR (2013) Bacteriophage 
genes that inactivate the CRISPR/Cas bacterial immune system. Nature 493: 
429-^32. doi:10.1038/naturell723. 



32. Seed KD, Lazimki DW, Calderwood SB, Camilli A (2013) A bacteriophage 
encodes its own CRISPR/Cas adaptive response to evade host innate immunity. 
Nature 494: 489-491. doi:10.1038/natarel 1927. 

33. Sheludchenko MS, Huygens F, Hargreaves MH (2010) Highly discriminatory 
single-nucleotide polymorphism interrogation of Escherichia coli by use of allele- 
specific real-time PCR and eBURST analysis. Appl Environ Microbiol 76: 
4337-4345. doi: 10. 1 1 28/AEM.OO 1 28- 10. 



PLOS ONE I www.plosone.org 



7 



June 2014 I Volume 9 | Issue 6 | e98811 



